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[57] ABSTRACT 

The physical realization of new solutions of wave prop¬ 
agation equations, such as Maxwell’s equations and the 
scaler wave equation, produces localized pulses of wave 
energy such as electromagnetic or acoustic energy 
which propagate over long distances without diver¬ 
gence. The pulses are produced by driving each ele¬ 
ment of an array of radiating sources with a particular 
drive function so that the resultant localized packet of 
energy closely approximates the exact solutions and 
behaves the same. 

20 Claims, 11 Drawing Sheets 
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There have been several previous attempts to achieve 
ELECTROMAGNETIC OR OTHER DIRECTED localized transmission beyond this Rayleigh distance 
ENERGY PULSE LAUNCHER with conventional systems. The best known of these are 


The U.S. Government has rights to this invention 5 
pursuant to Contract No. W-7405-ENG-48 between the 
U.S. Department of Energy and the University of Cali¬ 
fornia, for the operation of Lawrence Livermore Na¬ 
tional Laboratory. 
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BACKGROUND OF THE INVENTION 
The invention relates generally to transmission of 
pulses of energy, and more particularly to the propaga¬ 
tion of localized pulses of electromagnetic or acoustic 
energy over long distances without divergence. 15 

As the Klingon battle cruiser attacks the Starship 
Enterprise, Captain Kirk commands “Fire photon tor¬ 
pedoes”. Two darts or blobs of light speed toward their 
target to destory the enemy spaceship. Stardate 1989, 
Star Trek reruns, or 3189, somewhere in intergalactic 20 
space. Fantasy or reality. The ability to launch localized 
packets of light or other energy which do not diverge as 
they travel great distances through space may incredi¬ 
bly be at hand. 25 

Following the pioneering work of J. N. Brittingham, 
various groups have been actively pursuing the possibil¬ 
ity that solutions to the wave equation can be found that 
allow the transmission of localized, slowly decaying 
pulses of energy, variously described as electromag- 3Q 
netic missiles or bullets, Bessel beams, transient beam 
fields, and splash pulses. These efforts have in common 
the space-time nature of the solutions being investigated 
and their potential launching mechanisms, pulse-driven 
antennas. 35 

Brittingham’s original work involved a search, over a 
period of about 15 years, for packet-like solutions of 
Maxwell’s equations (the equations that describe how 
electromagnetic waves propagate). The solutions 
sought were to be continuous and nonsingular (well- 40 
behaved, realizable), three-dimensional in pulse struc¬ 
ture (localized), and nondispersive for all time (faith¬ 
fully maintaining their shape). They were also to move 
at the velocity of light in straight lines and carry finite 
electromagnetic energy. The solutions discovered, 45 
termed focus wave modes (FWMs), had all the afore¬ 
mentioned properties except the last; like plane-wave 
solutions to the same equations, they were found to 
have finite energy density but infinite energy, despite all 
attempts to remove this deficiency, and thus are not 50 
physically realizable. 

Conventional methods for propagation of energy 
pulses are based on simple solutions to Maxwell’s equa¬ 
tions and the wave equation. Spherical or planar wave¬ 
forms are utilized. Beams of energy will spread as they 55 
propagate as a result of diffraction effects. For a source 
of diameter D and wavelength of \ the distance to 
which a pulse will propagate without substantial spread 
is the Rayleigh length D 2 /\. 

Present arrays are based on phasing a plurality of 60 
elements, all at the same frequency, to tailor the beam 
using interference effects. In a conventional antenna 
system, such as a phased array driven with a monochro¬ 
matic signal, only spatial phasing is possible. The result¬ 
ing diffraction-limited signal pulse begins to spread and 65 
decay when it reaches the Rayleigh length L r. For an 
axisymmetric geometry, an array of radius a, and a 
driving wavelength of \, L r is about a 2 /\. 


the super-gain or super-directive antennas, where the 
goal was to produce a field whose amplitude decays as 
one over the distance from the antenna, but whose an¬ 
gular spread can be as narrow as desired. There are 
theoretical solutions to this problem, but they turn out 
to be impractical; the smallest deviation from the exact 
solution completely ruins the desired characteristics. 

The original FWMs can be related to exact solutions 
of the three-dimensional scalar wave equation in a ho¬ 
mogeneous, isotropic medium (one that has the same 
properties at any distance in all directions). This equa¬ 
tion has solutions that describe, for example, the famil¬ 
iar spherical acoustic waves emanating from a sound 
source in air. 

The FWMs are related to solutions that represent 
Gaussian beams propagating with only local deforma¬ 
tion, i.e., a Gaussian-shaped packet that propagates with 
changes only within the packet. Such a pulse, moving 
along the z axis, with transverse distance denoted by p, 


**('>')= e<k(z+c,) [ 

is an exact solution of the scalar wave equation devel¬ 
oped by applicant. This fundamental pulse is a Gaussian 
beam that translates through space-time with only local 
variations. These pulses can also form components of 
solutions to Maxwell’s equations. 

These fundamental Gaussian pulses have a number of 
interesting characteristics. They appear as either a 
transverse plane wave or a particle, depending on 
whether k is small or large. Moreover, for all k they 
share with plane waves the property of having finite 
energy density but infinite total energy. 

Thus traditional solutions to the wave equation and 
Maxwell’s equations do not provide a means for launch¬ 
ing pulses from broadband sources which can travel 
desirable distances without divergence problems. The 
laser is a narrowband light source which has a relatively 
low divergence over certain distances (i.e. relatively 
long Rayleigh length). However, acoustic and micro- 
wave sources, because of longer wavelengths, are more 
severely limited. Phased arrays do not provide the solu¬ 
tion. 

Accordingly, it is an object of the invention to pro¬ 
vide method and apparatus for launching electromag¬ 
netic and acoustic pulses which can travel distances 
much larger than the Rayleigh length without diver¬ 
gence. 

It is also an object of the invention to provide method 
and apparatus for launching pulses which approximate 
new solutions to the scalar wave and Maxwell’s equa¬ 
tions. 

It is another object of the invention to physically 
realize new solutions to the scalar wave and Maxwell’s 
equations which provide localized packets of energy 
which transverse large distances without divergence. 

It is a further object of the invention to provide com¬ 
pact arrays for launching these pulses. 
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SUMMARY OF THE INVENTION 
The invention is method and apparatus for launching 
electromagnetic and acoustic energy pulses which 
propagate long distances without substantial diver- 5 
gence. A preferred embodiment of the invention is 
based on the recognition that a superposition of the 
FWM pulses can produce finite-energy solutions to the 
wave equation and to Maxwell’s equations. As with 
plane waves, the infinite-energy property is not an in- 10 
surmountable drawback per se. The variable k in the 
solution provides an added degree of freedom, and 
these fundamental Gaussian pulse fields can be used as 
basis functions, a superposition of which represent new 
transient solutions of the wave -equation. In other 55 
words, these infinite-energy solutions can be added 
together, with the proper weighting, to yield physically 
realizable, finite-energy solutions. More generally, the 
invention applies to any nonseparable space-time solu¬ 
tion <t>*(r,t) of the relevant wave propagation equation, 2 o 
and may in some cases even be based on an approximate 
solution. 

For example, either the real or imaginary part of the 
function 

25 

Ar.t) = /“ <t>k(r,t)F(k)dk 


(where <I>* is the exact solution or an approximation 
thereof) is also an exact, source-free solution of the 3 
wave equation. The F(k) function is the weighting func¬ 
tion (the spectrum), and the resulting pulses having 
finite energy if F(k) satisfies certain integrability condi¬ 
tions. This representation utilizes basis functions that 
are localized in space and, by their very nature, are a 35 
natural basis for synthesizing pulse solutions that can be 
tailored to give directed wave energy transfer in space. 

A bidirectional representation is also possible, which 
leads to analogous solutions in geometries that have 
boundaries (propagation of waves in waveguides). 40 
Solutions to Maxwell’s equations follow naturally 
from these scalar wave equation solutions. Such elec¬ 
tromagnetic pulses, characterized by their high direc¬ 
tionality and slow energy decay, are called electromag¬ 
netic directed-energy pulse trains (EDEPTs). They are 45 
a step closer to a classical description of a photon, a 
finite-energy solution of Maxwell’s equations that ex¬ 
hibits a wave/particle duality. The corresponding 
acoustic pulses, which are solutions of the scalar wave 
equation, are called ADEPTS. 50 

In most general terms, the invention starts with a 
solution 4>jt( ,t) of the relevant wave equation (scaler 
wave equation, Maxwell’s equations or other equa¬ 
tions), chooses an appropriate weighing or spectrum 
function F(k), computes a drive function f(r,t) and ap- 55 
plies the appropriate drive function to each element of 
an array to launch pulses of wave energy. 

The invention particularly applies to broadband 
sources such as acoustic and microwave sources. Each 
element of an array of radiating elements is driven by 60 
the appropriate driving function for that individual 
element. The array is preferably a finite planar array, 
and may be folded to produce a more compact configu¬ 
ration. 

65 

BRIEF DESCRIPTION OF DRAWINGS 
FIGS. 1A and B show surface plots and correspond¬ 
ing contours of a Gaussian electromagnetic pulse with 


4 

Zo=l cm and k=0.333/cm at the initial position t=0, 
pulse center at z=0, and at t=7r ms with pulse center at 
z=942 km, respectively. 

FIGS. 2A and B show the energy density of a spread 
Gaussian pulse (small k, wavelike) and localized Gauss¬ 
ian pulse (large k, particle-like), respectively. 

FIGS. 3 A and B show the field energy density of the 
electromagnetic MPS pulse with a= 1 m, a= 1, 
b = 10 14 /m, /3=6X 10 15 , and z<,= 1 cm, normalized to its 
maximum value at t =0, at t =0, z =0 and at t= K^sec, 
z=9.42x 10 9 km, respectively. 

FIGS. 4A and B compare a traditional Gaussian 
beam solution of the wave equation with a localized 
transmission solution of the invention which has a large 
bandwidth and maintains its spatial and temporal fre¬ 
quency distribution over long distances. 

FIG. 5 is one quadrant of a schematic ADEPT- 
/EDEPT planar array of radiating sources, each driven 
with a particular time function having a large frequency 
bandwidth, with three representative driving pulses, 
each different, shown for three individual radiating 
elements. 

FIG. 6 illustrates how pulses radiated by each ele¬ 
ment in the array combine to form a resulting localized 
packet of wave energy. 

FIGS. 7A, B illustrate a drive function applied to the 
center of an array, and its frequency spectrum. 

FIGS. 8A, B illustrate a drive function applied to a 
non-center element of an array, and its frequency spec¬ 
trum. 

FIGS. 9A, B illustrate a more complex drive func¬ 
tion, and its frequency spectrum for a more noncentral 
element. 

FIG. 10 is a plot of field value or a function of array 
radius for several axial distances. 

FIG. 11 A, B illustrate the mapping of a large array 
into a folded array. 

FIG. 12 illustrates the ratio of the reconstructed to 
exact field value along the direction of propagation for 
20,000 element staggered and unstaggered 1.0 m arrays. 

FIGS. 13 A and B show the center drive function, 
and its frequency spectrum, for an experimental 
ADEPT array. 

FIG. 14 shows an experimental system to verify lo¬ 
calized transmission of ADEPT waves. 

DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENT 

The invention is method and apparatus for launching 
localized pulses of energy which substantially approxi¬ 
mate EDEPTs, electromagnetic directed energy pulse 
trains, which are exact pulse solutions of Maxwell’s 
equations in an isotropic, homogeneous medium, or 
ADEPTs, acoustic directed energy pulse trains, which 
are exact pulse solutions of the acoustic (scalar) wave 
equation in an isotropic, homogeneous medium. In a 
preferred embodiment, both classes of solutions can be 
constructed from “Focus Wave Modes”, which are 
exact solutions that represent Gaussian beams translat¬ 
ing through space with only local deformations. How¬ 
ever, the principles of the invention apply to any wave¬ 
form 0*(r,t) which is a nonseparable space-time solu¬ 
tion to the relevant wave equations. The relevant wave 
propagation equations are typically the scalar wave 
equation and/or Maxwell’s equations, but in some cases 
other wave equations may apply. Approximate solu- 
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tions may in some cases be substituted for exact solu¬ 
tions. 

If c is the propagation speed of waves in this medium, 
the axisymmetric solution (p 2 =x 2 +y 2 ) 


* kf.nl) = 


e -k p 2/[u+i(z-ct)] 
4iri[z 0 + i(z — c/)] 


(1) 


is a moving, Gaussian pulse modulated by a plane wave. 10 
It can have either a plane wave or a particle-like charac¬ 
ter depending on whether k is small or large. An added 
degree of freedom has been introduced into the solution 
through the variable k, and these fundamental Gaussian 
pulse fields can be used as basis functions to represent 
new transient solutions of the wave equation and to 
launch pulses. In particular, if s(p,z,t)=p 2 /[zo+i(z—ct)- 
]—i(z+ct), the function 

20 

-<*= (2) 
Ant) = J <S>lfr,l)F{k)dk = 

0 


+ ’(;-«)] /“ we -^>. 25 

is also an exact source-free solution of the wave equa¬ 
tion and provides the drive functions for a source array. 

In contrast to plane wave decompositions, it uses basis 3 
functions that are localized in space-time and hence, by 
their very nature, are a natural basis for synthesizing 
finite energy pulse solutions that can be tailored to give 
directed wave energy transfer in space. In particular, 35 
the representation yields finite energy solutions of the 
scalar wave equation if F(k) satisfies 


/" dk\F(k)\ 2 e lks <>Ei(2k2 0 ) < -±- /" dk\m\ 2 J~ < 00 ^ 40 

where E„(x) is the exponential integral function of order 
n. This occurs, for example, if k - lF(k) is square integra¬ 
te. 45 

FIGS. 1A and B show scaled plots of the fundamen¬ 
tal Gaussian pulse of Equation (1). The spatial variation 
about the pulse center at t=0, z=0 and t=7rms, z=942 
km is shown; the pulse is moving in the positive z direc- 50 
tion at the speed of light. FIGS. 2 A and B show the 
energy density of the fundamental Gaussian pulse (real 
part of [47ri<J>i]2), respectively, for small k where the 
pulse looks like a transverse plane wave, and for large k 
where the pulse is very localized and looks like a parti- 55 
cle. 

A preferred embodiment of the invention for produc¬ 
ing pulses which closely approximate the theoretical 
solutions is based on the spectral function ^ 

F{k)=4m{3(J3k-b) a - l e -«<0*-*>H(k-.b//8)/T(a) 

derived by scaling and truncating a power spectrum, 
where H(x) is Heaviside’s function and T(x) is the g 5 
Gamma function. This spectrum yields the associated 
modified power spectrum (MPS) pulse (drive function): 


Ant) = Re 



(4) 


where a and /3 are parameters that are chosen to 
achieve specific characteristics, and s is defined as be¬ 
fore. 

The physical characteristics of the MPS pulse are 
very appealing. This pulse can be optimized so that it is 
localized and its original amplitude is recovered out to 
extremely large distances from its initial location. In 
particular, for a distance z < </3/2b and z </8a/2, the 
amplitude of the pulse at the pulse center is constant. It 
then becomes oscillatory with an oscillation length of 
ir/b in an intermediate zone, /3/2b <z </Ja/2, recover¬ 
ing its original amplitude when z=n (7r/3/b), n being 
any positive integer. Finally, when the observation 
point is very far away from the origin, z>/3a/2, the 
MPS pulse decays like l/z°. Therefore, the initial am¬ 
plitude of the MPS pulse is recovered until the distance 
z~/Ja/2., and since j3 is a free parameter, this distance 
can be made arbitrarily large. The transverse behavior 
of this MPS pulse at the pulse center is essentially f(p,z 
=ct)~exp(—bp 2 //3z) f(0,z=ct). Thus, by adjusting 
the ratio b/j3z 0 , one can adjust the degree of transverse 
localization. The MPS pulse is also localized longitudi¬ 
nally, decaying along z as l/[ 0 2 +(z—ct) 2 ]away from 
the pulse center. 

To produce localized transmission (high directional¬ 
ity) and slow energy decay, the parameters a, b, o, y3B 
, and z<,of the MPS pulse are selected to achieve a pulse 
within the microwave spectrum and possibly within the 
realm of our physical appreciation and experience. In 
particular, set the MPS parameters to be a=lm,a=l, 
b= 10 14 /m, y3=6x 10 15 , and z 0 = 1 cm. The peak of the 
spectrum of this pulse is in the microwave region at 8.4 
GHz. (However, these parameters can be varied to 
design pulses with similar characteristics in different 
frequency regimes.) 

FIGS. 3A and B show surface plots and the corre¬ 
sponding contour plots of the electromagnetic energy 
density of the electromagnetic MPS pulse relative to the 
pulse center locations z=0 km and 9.42X 10 9 km (9.42 
billion kilometers, farther than from Earth to Pluto). 
These results definitively show the localization of the 
field near the direction of propagation over very large 
distances. The recovery of the initial energy density in 
this case actually occurs out to z=/8a/2=3.00X 10 12 km 
(about one-third light year). 

A comparison of this specific choice of pulse repre¬ 
sentation and the classically popular Gaussian beam 
shows that the electromagnetic MPS pulse characteris¬ 
tics are dramatically better. As shown in FIGS. 4 A and 
B the waist of zero-order Gaussian beam field at z=0 is 
w and X is its wavelength. Along the direction of propa¬ 
gation its amplitude varies as l/[l-f(\z trw 2 ) 2 ]! so that 
the distance to the near/far-field boundary or Rayleigh 
length, where it begins to decay as 1/z, is nominally 
reached when z~?rw 2 A. The square of the amplitude 
at that point is half its initial value and its radius has 
spread to VT" w and increases (diverges) as 
9z ~ (X/7rw)z in the far-field. 

With the above MPS pulse parameters, the waist w of 
an equivalent Gaussian beam=(b//3z 0 )-1=1.3 m. The 
highest frequency of the spectrum for this MPS pulse is 
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f max =50 GHz; the corresponding wavelength is 
\ m i„=c/f m ax =6 nun. For the equivalent diffraction- 
limited Gaussian beam field these defining parameters 
give the distance to the far-field as 0.872 kin and the 
spread of the field at 10 10 km as 1.5 X10 7 km. The field 5 
amplitude at 10>°km is essentially 10— 10 its initial value. 
The localization of the electromagnetic MPS pulse near 
the z-axis and the recovery of its initial amplitude well 
beyond the classical far-field distance confirms that the 
MPS pulse has propagation characteristics that are 
much better than the corresponding diffraction-limited 
Hermite-Gaussian laser field. 

The generation of pulses of energy that satisfy the 
wave equation and provide localized transmission of 
energy requires an antenna, or finite planar array of 
radiating elements, which produces fields that are sub¬ 
stantially similar to the exact mathematical solutions so 
that the resulting pulses behave substantially as pre¬ 
dicted. 

Ideally, the antenna system is a finite planar array of 
point sources each of which radiates spherical pulses 
that can be combined using a Huygens representation 
into the array field. (The Rayleigh distance for a point 
source is zero, so we are always in the far field of each 
radiating element and can readily obtain the overall 
field response of the array by superposition.) 

The antennas, as illustrated in FIG. 5, include circu¬ 
lar, rectangular, and hexagonal arrays of equally spaced 
elements. The arrays may be planar or nonplanar. The 
driving function for each element is a broad-bandwidth 
waveshape determined from the exact wave-equation 
solution and its derivatives. This is marked contrast to 
conventional arrays, whose elements are driven with 
monochromatic signals. The elements include ultra¬ 
sonic acoustic transducers (piezoelectric) and micro- 
wave sources such as dipoles, horns, etc. The invention 
can also be applied to lasers and other sources. 

As shown in FIG. 6, the resulting field (the sum of 
these individually radiated time histories) is a localized 
pulse that maintains its shape and compactness at dis¬ 
tances well beyond the conventional Rayleigh distance. 
Furthermore, the ADEPT/EDEPT-driven arrays ap¬ 
pear to be very robust (not strongly sensitive either to 
parameters defining the array or to perturbations in the 
initial aperture distributions). 

The elements in the array are driven by a waveform 
determined by Eq. (4) for the particular case of the MPS 
pulse. The drive functions vary with the positions of the 
elements in the array. In the more general case, for any 
selected spectral function F(k), the drive function for 
each element in the array is determined by Eq. (2). A 
digital waveform synthesizer can be utilized to generate 
the appropriate waveforms. FIGS. 7A, B illustrate the 
drive function at the center of an EDEPT array, and its 
Fourier spectrum; the spectrum shows the broadbound 
nature of the excitation. FIGS. 8A, B illustrate the drive 
function and spectrum applied to another source not at 
the center; FIGS. 9A and B are a more complex wave¬ 
form, and its spectrum, applied to a more noncentral 
source. The spectrum of FIG. 7B forms the envelope of 
the spectra for all the other driving functions; e.g., the 
spectrum of FIG. 8B fits within the spectrum of FIG. 
7B as does the spectrum of FIG. 9B. 

The larger the array, the closer the physically realize- 
able pulse approximates the theoretical solution. To 
produce MPS pulses with a rectangular array of (N+1) 
X(M+1) equally spaced elements, the drive functions 
are approximated by 


Ar.t) = - 

W n * x ’ miky -*’‘ ~ R nm/c)*xAy] A J Rnm . 

where the array element spacings are Ax and Ay and the 
distances 

Rrnn=[(x-n &x) 2 +(y-m Aji) 2 +(z-/) 2 ]i 

For a finite circular array with distinct, equally spaced 
annular sections, the drive function is approximated by 


Ap = <U/) = - - R n / c)Area(n)] . 

where the distances p' n = n A p' and 

R«=(i>>' n 2 -Kz-z')2]i 

and the area weighting Area (n)=2 irp„' A p'. The func¬ 
tion 'P is related to the function 0 and is calculated from 
the drive function f by 

af/az , -[af/acf](2-/)/it-A2-z')/« 2 - 

Nonuniform spaced arrays can also be used. 

A quantitative measure of the size of a circular array 
needed to reconstruct the MPS pulse at increasingly 
larger axial distances is shown in FIG. 10. The field 
values at the axial distances z=l km, 10 km, 100 km, 
1000 km and 10,000 km generated by a finite circular 
array are plotted against the array radius. The curves 
become horizontal when the exact field value is 
reached. The size of the array controls the reconstruc¬ 
tion distance. The scaling ratio is approximately a factor 
of 10 increase in radius size (in cm) for a factor of 10 
increase in distance along the direction of preparation 
(in km). A smaller array does not increase the transverse 
width of the pulse, but only degrades the reconstruc¬ 
tion. As the array size is increased, the pulse definition 
is enhanced relative to the surrounding fields. 

A difficult issue, and probably the one most used as a 
figure of merit is the distance over which localization 
will occur for the ADEPT/EDEPT array-launched 
pulse as compared with traditional fields. There is no 
exact value for the Rayleigh distance in the case of an 
ADEPT/EDEPT array because of the broadband 
pulsed nature of these solutions, which have little in 
common with the monochromatic radiation on which 
the Rayleigh-distance concept is based. Whether the 
Rayleigh distance is dramatically surpassed or is simply 
reached by the new array may be moot—it depends 
greatly on the intended application. However, even 
from a modest size array, energy can be transmitted 
locally without spreading over significantly large dis¬ 
tances. 

Array reconstructed MPS pulses are not very sensi¬ 
tive to perturbations in the initial driving functions. An 
amplitude taper, e.g. a Hanning window, applied to the 
aperture driving functions actually helps the pulse re¬ 
construction by decreasing the late time oscillations and 
the source density, with a decrease in peak amplitude. 
With a slightly larger aperture size and the amplitude 
taper, a much improved reconstruction of the MPS 
pulse is achieved. The effect of the taper is to remove all 
of the low frequency components of the field and to 
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slightly emphasize the high frequency ones. The effect 
of frequency filtering the driving functions is also posi¬ 
tive. Late-time oscillations occur because of the pres¬ 
ence of undesirable, ill-behaved higher frequency com¬ 
ponents which become emphasized in the superposition. 
The use of a filter to remove some of these unwanted 
components of the initial driving functions, e.g., a low 
pass 2nd order Butterworth filter, removes much of the 
late time noise at a slight cost in the peak value. The 
application of random Gaussian noise to the individual 
driving functions also has little adverse effect. 

One technique of obtaining a large array reconstruct 
from a smaller array is to use a folded array, as shown 
in FIGS. 11A, B. The exterior of a planar circular array 
of radius t max is folded onto its interior with the confor¬ 
mal map p—*r ma x 2 /p- The folded array may be stag¬ 
gered or unstaggered. The folded array is staggered if, 
when the interior points are located at nAp, the mapped 
exterior points are at (n+£)Ap; unstaggered if the inte¬ 
rior and mapped exterior points coincide. The folding 
trades a less complicated source distribution for a more 
complex one. The drive functions for all points in the 
folded array are determined by their positions in the 
unfolded array. As an example, the complex drive func¬ 
tion of FIGS. 9A, B would be applied to an element 
which has been folded in near the center of the array. 

The ratio of the reconstructed to exact field valve 
along the direction of preparation for 20,000 element 
staggered and unstaggered l.Om arrays are shown in 
FIG. 12. A conventional single frequency phased array 
would only reach to about 10 5 m at best, whereas the 
arrays of the invention reach beyond 10 8 m, an improve¬ 
ment by a factor of at least 1000. 

The invention has been demonstrated in the acoustic 
regime based on ultrasonic pulses propagating in water. 
Acoustic directed energy pulse trains (ADEPTs) are 
localized wave packets which are based on novel solu¬ 
tions to the scaler wave equation which are unique in 
their intrinsic space-time nature. 

The ADEPTs can be generated with a finite array of 
radiating elements by specifying both their spatial and 
their temporal distributions. The driving functions for 
the array elements are determined by the exact solution 
and its derivatives. Computer simulations, in particular, 45 
a Huygens reconstruction based on a computer model 
of a finite planar array of point sources reproduces the 
modified power spectrum (MPS) pulses at large dis¬ 
tances away from the array. The array-generated MPS 
pulse appears to be very robust and insensitive to per- 50 
turbations in the specified source distributions. These 
results are also insensitive to the type of array (circular, 
rectangular, or hexagonal) considered. The feasibility of 
launching an acoustic MPS pulse travelling at the speed 
of sound in water (1.5 km/s) is demonstrated expert- 55 
mentally. The ADEPT parameters selected are a=lm, 
b=600 m, a= 1,0=300 and z o =0.45 mm. FIGS. 13 A, 

B show the associated driving function, and its fre¬ 
quency spectrum, which is applied to the center ele¬ 
ment of the array. 

The pulse parameters are chosen to facilitate the 
design of an acoustic experiment. Its pulse (waist =1.5 
cm) and its spectrum (peak at 0.6 MHz and practically 
all of its energy below 2.0 MHz), coupled with the 
choice of array geometry and array element size, permit 65 
the desired effects to occur in a distance less than 2 m, 
the effective length of an available water tank. In partic¬ 
ular, with the Rayleigh distance Lj? at the peak fre- 


10 

quency being 28 cm, there are about 7L r available in 
the tank. 

The pulses are generated with an array that is compu¬ 
tationally and experimentally simple and within the 
5 scope of limited experimental resources. The latter limi¬ 
tation imposed some significant constraints and dictated 
the experimental arrangements. 

To produce the field generated by an array, it would 
be necessary to drive each element of the array with the 
10 appropriate waveform for that particular element. Thus 
a generator and transducer for each element would be 
needed. 

In the first stage of the experimental program, the 
complexity introduced by independent generators and 
15 transducers is minimized by constructing a synthetic 
array. In such a construct, each element of the array is 
driven individually by its own source, each source hav¬ 
ing an appropriate waveform; the field generated in this 
configuration is recorded. This procedure requires only 
20 one generator and one transducer (element) at any one 
time to provide the contribution of a given element to 
the array performance. After all elements are so driven 
and their radiated fields recorded, the array field is 
synthesized by superposition of the fields previously 
25 recorded. Thus, it is possible to generate the field radi¬ 
ated by an array by computational reconstruction, using 
experimentally measured contributions from individual 
radiators. 

An acoustic field detector was used that has a mini- 
30 mal impact on the measured field quantity, i.e., the 
presence of the detector does not impact the variable 
being measured. A laser beam/photodiode combination 
that measures sound-wave induced changes in the re¬ 
fractive index of water was used. The laser beam propa- 
35 gates through the water tank at right angles to the direc¬ 
tion of sound propagation and through a window to 
strike a photodiode. A lens placed between the water 
tank and the photodiode puts a reduced virtual image of 
the photodiode within the tank near the sound field. 
40 The sound field creates local variations in the optical 
refractive index of the water in the tank. At low ultra¬ 
sonic frequencies, the total effect on the optical beam is 
a phase modulation that occurs at a unique plane as the 
light beam transverses the sound wave. 

The magnitude of the phase modulation is the effec¬ 
tive optical path length through the local variations of 
refractive index. It is the line integral of the refractive 
index through the regions of optical retardation and 
speed-up induced by the sound. The intensity of the 
light propagating across the sound field is proportional 
to the curvature in the wavefronts, which is, in turn, the 
second derivative of the phase modulation. The photo- 
detector thus measures the second time derivative of the 
acoustic field interacting with the laser beam. 

Primary interests have been directed toward measur¬ 
ing the fields radiated by an array at a point in space, not 
along a line, but the detector system measures line inte¬ 
grals of the field. In order to circumvent this problem, 
the reciprocity principle is used, by which the role of 
60 transmitter and receiver can be interchanged without 
affecting the measured response. That is, the principle 
of reciprocity is used to interpret the field radiated by a 
point source and measured by a line detector as the field 
radiated by a line source and measured by a point detec¬ 
tor. 

In summary the array is one-dimensional and syn¬ 
thetic and uses reciprocity to reverse the roles of sender 
and receiver. The practical result is that one can deter- 
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mental Gaussian control and the corresponding numeri¬ 
cal simulations, a Gaussian beam field was fit to the 
data. The effective frequency was 0.6 MHz, the peak of 
the spectra of the driving wavefunctions. The effective 


mine the wave that each element of a pulsed array of 
line sources must radiate to produce an ADEPT pulse 
by substituting line detectors and irradiating them with 

waves emanating from point sources that mimic the - r - 

properties of the ADEPT wave. By synthesizing the 5 Rayleigh length for the experiment was thus about 28 
array, one can use just one detector and move it to 
different positions. For each position of the detector a 
time history of the received pulse is recorded. Then all 
the different time histories are added up to yield the 


The synthetic linear array experiment is simulated 
(theoretical ADEPT) by driving each element at (x,y) 
in a rectangular array with the wavefunction at (O.y). 


radiation the array would emit if we were pulsing its 10 This ensures that the array appears “linear,” as does the 


elements, each with the wave shape received at its loca¬ 
tion in the array. 

The experiment system used to verify localized trans- 
i of ADEPT waves is shown in FIG. 14. The 


reciprocal laser diagnostic system. The simulated array 
is 6 cmX6 cm and contains 441 elements in a 21x21 
equally spaced pattern. 

Experiments and computer simulations verify the 
sending system consists of a progammable waveform 15 existence, behavior and practical realizability of the 


generator (pulser) 20, a power amplifier 22, and a piezo¬ 
electric ultrasonic transducer 24, with a command com¬ 
puter 26 to down-load the different waveforms to the 
generator 20. The receiving system consists of a laser 


ADEPT solutions with the parameters described 
above. The experimental and theoretical ADEPT re¬ 
sults are in excellent agreement. Both show the remark¬ 
able retention of the pulse behavior well beyond two 


beam 28 from CW HeNe laser 30 to probe the ultra- 20 R ay j e i g h distances. The contrast in compactness be¬ 


tween the ADEPT pulse and a Gaussian pulse is also 
evident. Both experimentally and numerically the linear 
array produces fields that begin to break up after 50 cm, 
about 2 Ljj. The rectangular array generated ADEPT 


sound pulse field 32 and a lens 34 -photodiode 36 combi¬ 
nation that places a virtual image 38 of the photodiode 
inside the immersion tank 40 very near the critical 

plane. The position controller 42 moves the transducer auuul A __ 

24 forward and backward, up and down, and sideways 25 avoids this effect because the off-axis elements compen- 
to place the virtual image 38 m different parts of the satg ginct splitting, 
pulse field 32. The system computer 44 issues com- - ■ - 

mands to the pulser 20 through command computer 26 
and to the position controller 42 and files data from the 


Even a simple unoptimized array can reproduce the 
MPS pulse at significant distances. More optimum de¬ 
signs will reach far beyond the classical Rayleigh 


photodiode 36 received through digitizer 46. Sync- 30 ^ advanced design squeezes a large 

nnfser 48is connected to nulser 20and through a delav ® . .. ~ , • .... _ 


pulser 48 is connected to pulser 20 and through a delay 
to digitizer 46. 

The sound source consists of a single, commercial 
ultrasonic transducer designed for nondestructive test¬ 
ing. The transducer is a piezoelectric disk 6.2 mm in 35 
diameter, with acoustically matched damping material 
on its backside, that produces a piston-like motion in the 
water. For distances greater than 6 cm, the resulting 
sound beam is in the far field of the transducer. Within 
the approximations of the design, it is a signal propor- 40 
tional to the third derivative (one from the transducer 
and two from the optics) of the driving function (the 
electrical signal applied to the source element) that is 
eventually acquired by the data acquisition system. 


array into a smaller one. One way is to drive the array 
with a more complicated set of pulse shapes whose 
functions are derived from the ADEPT solution within 
the array distance and beyond it. The solutions beyond 
the ma ximum distance are “folded” into the interior of 
the array in a particular way, trading a simple source 
distribution for a much more complicated one. 

The results of this maneuver on the array-launched 
ADEPT field is dramatic. Localization for the above 
case out to about 10 m is obtained which is more than 30 
L R. 

The unique aspect of the ADEPT/EDEPT solutions 
is their intrinsic space-time nature. An MPS pulse c 


In the coordinate system indicated in the 45 be designed to recover it initial amplitude after propa- 

- - - very large distances while spreading very little. 


linear array of the experimental element positions 
along the y-axis, the laser measurement is taken along 
the x-axis, and the direction of propagation is along the 
z-axis. The synthetic linear array consisted of 21 ele- 


The pulse moves virtually unchanged in the “near” 
zone, “sloshes” about the pulse center in the “intermedi¬ 
ate” zone, recovering its initial amplitude at intervals 


positions symmetrically arranged about y=0 with 50 out to very large distances, and finally falls off as in- 


3-mm separation. Thus, 11 driving functions 
ployed, and the total array width was 6 cm. Time 
cords of length 10.24 fis were taken. The 176 unique 
experimental waveforms, launched with a normalized- 


verse distance in the “far” zone. These pulses can be 
produced with a finite array of radiating elements indi¬ 
vidually driven with appropriately shaped pulses. A 
Huygens reconstruction based on the causal, time- 


to-one maximum amplitude, were weighted and super- 55 retarded Green’s function and^ a^finite planar array^ of 
posed to construct the field of the synthetic linear array. ~ * r '” ,c " Urrro H,c - 

It should be emphasized that, save for multiplying the 
waveforms by a weighting coefficient dependent on the 
assumed position in the array for each generated signal, 
no other processing was applied, either before or after 60 
the experiment. 

As a control, a Gaussian beam was constructed from 
the same array with different weights but the same 
waveform from all array elements. Each element of that 


point sources reproduced the MPS pulses at large dis¬ 
tances. The array-generated MPS pulse appears to be 
very robust and insensitive to perturbations in the initial 
source distributions. 

The physical realization of new solutions of Max¬ 
well’s equations and the wave equation provide the 
possibility of propagating localized pulses of electro¬ 
magnetic or acoustic energy over long distances with¬ 
out loss. Such localized transmissions could have appli- 


array was driven with the center ADEPT driving func- 65 cations in communications, remote sensing, power 


tion and was weighted with a Gaussian amplitude 
exp(—p 2 /w 0 2 ), where w 0 = 1.5 cm, the same initial 
transverse waist as the ADEPT. For both the experi- 


transmission, and directed-energy weapons. 

Changes and modifications in the specifically de¬ 
scribed embodiments can be carried out without depart- 
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ing from the scope of the invention which is intended to 
be limited only by the scope of the appended claims. 

I claim: 

1. A method of producing a localized packet of wave 
energy which travels substantially large distances com- 4 5 6 7 8 9 
pared to the Rayleigh length without substantial diver¬ 
gence, comprising: 

independently driving each element of a finite array 
of radiating elements with a drive function deter- 10 
mined by 


14 

10. The method of claim 1 further comprising form¬ 
ing the drive functions with a digital waveform synthe¬ 
sizer. 

.11. Apparatus for producing a localized packet of 
wave energy which travels substantially large distances 
compared to the Rayleigh length without substantial 
divergence, comprising: 
a finite array of radiating elements; 
driving means associated with each radiating element 
of the array which apply a driving function 


fir.t) = 


r 


0 


QkHr.omdk 


15 


where <t>/t(r,t) is a basis function which is an exact or 
approximate non-separable space-time solution of 
the relevant wave propagation equation and F(k) is 
a spectrum function which satisfies 20 


Ar.t) = f°° *n{r.mk)dk 

o 

where 4>*(r,t) is an exact or approximate non-separable 
space-time solution of the relevant wave propagation 
equation and F(k) is a preselected spectrum function 
which satisfies 


S* dk \m\ 2 ^ he Ei(2kzo) < 
0 


-ST f dk \™?T 

2. The method of claim 1 wherein 


4>Ar,!) = 


,,, e -kp2/[ze+i(s-ct)) 

4iri[zo + <(z - «)] ’ 


f dk \m^ ke E X (2kz 0 ) < 
0 


25 


-W f dk Wf<“' 


12. Apparatus of claim 11 wherein 


30 


4>*(r,/) = e/k(s+ci) 


where s(p,z, t) = p 2 /[z 0 +i(z—ct)]—i(z+ct). 

3. The method of claim 2 wherein F(k) is a modified 
power spectrum 

F\k)=4ni0(fik-b) a -'i a ^k-b)H(k - b/p)/T(a) 

45 

where H(k—b//3) is Heaviside’s function and T (a) is 
the Gamma function, and 


35 **>«■ 4 wto + H-Sr f dkWK-^\ 

where s(p,z,t) 32 

2 /[z<>+i(z—ct)]—i(z+ct)- 

40 13. Apparatus of claim 12 wherein F(k) is a modified 

power spectrum 

F[k)—4nip(fik—b) a ~ l ef l ^k—b)H(k—b/P)/r(,a.) 

where H(k—b//3) is Heaviside’s function and T(a) is 
the Gamma function, and 



4. The method of claim 1 comprising forming the 
radiating elements of broadband sources. 

5. The method of claim 4 comprising forming the 
radiating elements of piezoelectric transducers. 

6 . The method of claim 4 comprising forming the 
radiating elements of microwave sources. 

7. The method of claim 1 further comprising folding 
the finite array of radiating elements into a compact 
folded array of preselected maximum dimensions. 

8 . The method of claim 1 further comprising forming 
the array in a planar geometry. 

9. The method of claim 8 comprising forming the 
array in a circular, square or hexagonal array. 


14. Apparatus of claim 11 wherein the radiating ele¬ 
ments are broadband elements. 

15. Apparatus of claim 14 wherein the radiating ele¬ 
ments are piezoelectric transducers. 

16. Apparatus of claim 14 wherein the radiating ele¬ 
ments are microwave sources. 

17. Apparatus of claim 11 wherein the array is a 
folded array. 

18. Apparatus of claim 11 wherein the array is a pla¬ 
nar array. 

19. Apparatus of claim 18 wherein the array is a cir¬ 
cular, square, or hexagonal array. 

20. Apparatus of claim 11 wherein the driving means 
includes a digital waveform synthesizer. 
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[57] ABSTRACT 

The physical realization of new solutions of wave prop¬ 
agation equations, such as Maxwell’s equations and the 
scaler wave equation, produces localized pulses of wave 
energy such as electromagnetic or acoustic energy 
which propagate over long distances without diver¬ 
gence. The pulses are produced by driving each ele¬ 
ment of an array of radiating sources with a particular 
drive function so that the resultant localized packet of 
energy closely approximates the exact solutions and 
behaves the same. 
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